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Abstract Particle creation has been considered as a possible justification for 
the accelerated expansion of the universe, obeying the second law of thermo¬ 
dynamics, together with the possible existence of Dark Energy. This paper 
introduces the possibility that the destruction of baryonic and/or dark matter 
particles also verifies the second law of thermodynamics thanks to a particle 
exchange with dark energy. General equations for the variation of the number 
of particles in accelerated universes have been obtained. Finally, a new model 
of the universe has been developed which predicts dark energy properties as 
well as particle exchange processes between dark energy and baryonic and/or 
dark matter. 
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1 Introduction 

It seems well proven that the universe is experiencing an acceleration process 
[1-3]. Physicists have elaborated various answers [4] to explain this accelerated 
state, among them, the following may be mentioned: the existence of a new 
kind of energy (called ’dark energy’ that we could consider to include Einstein’s 
famous cosmological constant) [5], matter creation [6-9], and a whole set of 
theories that modify or go beyond the General Theory of Relativity, such as 
F(R) [10]. 
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The problem in elucidating what is really happening resides in the cur¬ 
rent data and the precision with which they are obtained. Many theories (and 
within them a whole set of different models) fit the available data, thus caus¬ 
ing a high degree of degeneration in solving the question of the accelerating 
universe. 

One of the ways of restricting this high degree of degeneration is by lay¬ 
ing down the conditions (theoretical and experimental) the different models 
should fulfill. In the research we are presenting, we have obtained the general 
equation that - from a thermodynamically point of view [11,12] - the uni¬ 
verse should fulfill in adiabatic conditions, when treated as an open system 
with particle creation or destruction [13-22]. From the previous equation, we 
have deduced the general equation a FLRW universe (treated as an adiabatic 
perfect fluid with particle creation or destruction) [23, 24] must verify by in¬ 
troducing the necessary terms in the energy-momentum tensor considering, 
as well, the possible existence of Dark Energy. In the present research, at no 
point the condition of constant specific entropy has been imposed [6, 24-26], 
only the condition of adiabatic open systems. 

Once the general condition an adiabatic FLRW universe with particle cre¬ 
ation has been established and in order to delimit the solutions of the FLRW 
general equations, we have deduced the general restriction these universes must 
satisfy in order to show an accelerated expansion such as the one our universe 
is currently experiencing. 

The rest of the article studies firstly and briefly the thermodynamics of 
open adiabatic systems with particle creation, then expanded Friedmann equa¬ 
tions for FLRW universes are established, and the thermodynamic restriction 
that allows Friedmann equations to adapt to open systems with particle cre¬ 
ation is set. In the next section, the general restriction an accelerated adiabatic 
universe with particle creation must satisfy is obtained. Later, two particular 
cases are analyzed: Prigogine model, and a new model developed for the first 
time (as far as the author is aware) in this article. To finish, the final conclu¬ 
sions are established. 


2 Thermodynamics of adiabatic open systems with matter creation. 

The first law of thermodynamics which applies to open systems (also called 
Gibbs equation for open systems) can be obtained from the first law of ther¬ 
modynamics applicable to closed systems (Eq.l) 


dE = dQ — pdV 

Defining the energy density (p) and the density of a particle (n)as 


( 1 ) 
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Where N and V are respectively the number of particles and the volume of 
the system. Applying equation 2 on equation 1, considering N constant, and 
defining dq = dQ/N it becomes 

d{-) = dq - pd{-) (3) 

n n 

Equation 3 can be applied to open systems, if it is considered that N may 
not be constant. Defining h as the enthalpy density, and supposing that N may 
not be constant, Gibbs equation for open systems is obtained from equation 3 
(Eq.4) 


d{pV) = dQ — pdV +—d{nV) (4) 

n 

Considering adiabatic systems {dQ = 0) equation 4 becomes (Eq.5) 

d(pV) + pdV ——d(nV) = 0 (5) 

n 

Expanding equation 5 and dividing everything by Vdt, it becomes (Eq.6) 

p = —fi (6) 

n 

Note that in order to deduce equation 6 it has not been necessary to impose 
the condition of a constant specific entropy [6, 24-26], just the condition that 
the system is adiabatic (although, as it will be shown later, this condition 
implies that the specific entropy is constant). 

Equation 6 can be applied to an adiabatic universe model where there 
is a variation in the number of particles either by their creation or by their 
elimination. 

In order to establish the condition that must be fulfilled by an adiabatic 
ELRW universe with particle creation, it is preferable to express equation 6 in 
a different way. 

Defining a covolume V = ka^, where k is a constant and a is the scale 
factor of the ELRW metric, it is easy to obtain (Eq.7) 


n 

n 



( 7 ) 


Where H is the Hubbel parameter. Substituting equation 7 into equation 
6, one gets (Eq.8) 


p + 3H(p + p) -—(p + p) = 0 (8) 

The study of the entropy variation of the universe provides very useful in¬ 
formation. In principle, entropy variation can come from two sources: the heat 
variation of the system, or the variation in its number of particles. Expressing 
this in the form of an equation, and using the definition of entropy, we have 
(Eq.9) 
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TdS = dQ + T^dN (9) 

Where T is the system temperature and the quotient ^ is the specific 
entropy a. 

Assuming that the universe is adiabatic, equation 9 is reduced to (Eq.lO) 

dS = adN (10) 

From equation 10, it follows that: 

1. The elimination of particles in the universe is not possible if the second 
law of thermodynamics is to be met. 

Dividing equation 10 by dt, one obtains S = aN. According to the sec¬ 
ond law of thermodynamics S' > 0, and since cr is a positive quantity, 
consequently fV > 0. 

2. Specific entropy is constant. 

Calculating the total derivative of the expression that defines specific en¬ 
tropy, S = Na, one obtains dS = Nda + adN, and comparing this expres¬ 
sion with equation 10, it follows that da = 0, therefore, cr is constant. 


3 Expanded FLRW equations, thermodynamic conditions for the 
creation of matter. 


The universe is assumed to be an ideal fluid in which metric is given by the 
metric of the FLRW universe (Eq. 11). 


ds"^ = -c^dt^ + a^{ _ „ -I- r'^dd'^ + r^Sin^0d<P^) (11) 

1 — eN 

In order to obtain more general Friedmann equations, two terms are added 
to the energy momentum tensor A first term (called A) is added to the 
energy term Too, and a second term (called B) is added to the pressure terms 
Tii. After carrying out the usual operations, a new set of Friedmann equations 
is deduced (Eq.l2) 


2 SttG (T c? 
= —arO+—A- — 


3c2 


■A 


2ff + 3H^ = -^p-^e + Ac^-B ( 12 ) 

Where c is the speed of light, G is the universal gravitation constant, and 
A belongs to the possible term attributed to dark energy. 

From the set of equations 12 it is easy to obtain the equation that allows 
us to deduce the condition a universe with matter creation must meet (Eq.l3) 

p + 3H{p + p) + [H (3A + B) + A = 0 


(13) 
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In the deduction of equation 13 , we have considered the general case where 
the terms A, B, and A may vary with time. 

By comparing equation 13 with its corresponding equation according to 
the thermodynamic development (Eq.8), we obtain the general condition that 
terms A and B must satisfy in order to have a universe with matter creation 
(Eq.l4) 


TV 3c^ ■ c^yi 

-]^(P+P) = + B) + A+—] (14) 

If dark energy density is considered to be constant, or if its existence is 
simply not considered at all, equation 14 becomes (Eq.l5) 

-^(P + P) = ^mSA + B) + A] (15) 

According to equation 8, the variation of density with time is 

N 

p=-3iJ(p + p) + —(p + p) (16) 

or, using equation 14 

P = —3i ?[p + p) — [H (3A + B) + A-\ —(17) 

From equation 16, defining the quotient ^ as B, and using the definition 
of enthalpy density h, we get 

p = h{r- 3H) (18) 

Taking into account that - = r — 3H, equation 18 is finally expressed as 
equation 19, which matches the one obtained by other authors [23] 

11 

p=h- (19) 

n 


4 Accelerated expansion of the universe. 

Since the universe is currently experiencing an expansion phase, it is advisable 
to lay down the condition for the existence of this acceleration, which will also 
allow for further delimiting the possible cosmological models. 

This condition has been obtained starting from the calculation of the ex¬ 
pression given by equation 20 

2^ > 0 (20) 

Whose sign will depend only on the second time derivative of the scale 
factor, since a is always positive. 

Starting from Friedmann equations (Eq.l2), and using equations 8 and 14, 
the general condition the accelerated universe must meet is obtained (Eq.21) 
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—+ 3p) H ——A — {A + B) > 0 (21) 

Eliminating the dark energy term, the final expression becomes (Eq.22) 

—+ 3p) — {A + B) >0 (22) 

It follows that the transition points from an accelerated phase to a decel¬ 
erated one must verify (Eq.23) 

(p + 3p) = {A + B) (23) 


5 Particular cases. 

In this section two specific cases of singular importance are studied starting 
from previously developed equations. 


5.1 Prigogine model. 


This model studied by Prigogine [24] and other authors [21, 23], assumes that 
parameters A and A are null, so equations 14 and 21 become 



(24) 

3^2 (p + 3p) B >0 

(25) 


Since H has the same dimensions as the quotient , it follows from equation 
24 that the term B can be interpreted as a pressure multiplied by some suitable 
dimensional constant. The constant is easily deduced from Eriedmann’s second 
equation [12] 


B = 


SttG 
2 


(26) 


Where Pc is the creation pressure caused by the creation of particles. After 
solving equation 24 for B, and comparing it with equation 26, we have 


h N 
“ ~mN 


(27) 


Erom equation 27 follows that matter creation creates a negative pressure. 
After some simple transformations, equation 27 can be expressed in the same 
form as the one found in [23] (equation 28) or in the form given by Prigogine 
[24] (equation 29) 
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h dLnin). 


(28) 


Pc = -- 


h d{nV) 
dV 


(29) 


The condition that particle creation velocity must meet in so that the uni¬ 
verse experiments an accelerated expansion is obtained by combining equations 
25, 26, and 27 (equation 30) 


N > , 30 ) 

If we assume that pressure follows equation of state of the p = Wp type, 
equation 30 becomes equation 31 which coincides with the one deduced in [21] 
when W = 0 


N > 


1 + 5W 

+ w 


(31) 


5.2 Transmutation model T 

This model, described for the first time in this article (as far as the author 
knows), is defined by the following parameter values: A = B = 0, A 0. 

With the aforementioned values, equations 14 and 21 are now expressed as 


TV 3c^ c^A 

J^ip + p) - 3 

(32) 

SttG 2c^ 

+ 3 A 

(33) 


Before analyzing the different submodels that can be deduced from equa¬ 
tions 32 and 33, it is necessary to analyze again equations 8, 10, 12 and 13 in 
the case dark energy exists. 

The first Equation 12 shows that when A = B = 0 the energy content of 
the universe is divided in three large groups: the group defined by the density p 
attributable to baryonic and dark matter, the curvature energy (which depends 
on the value of e), and finally the dark energy variable (which depends on the 
value of A). 

Equation 8 can be applied only to baryonic and to dark matter as reflected 
when comparing equations 8 and 13 in order to obtain equation 14. From what 
has been stated before, it can be observed that in equation 32 the number of 
particles N and its temporal variation N relate only to baryonic and dark 
matter. It should be also noted that if the pressure p corresponds to baryonic 
and dark matter then it follows that the pressure of dark energy would be zero 
(which is consistent with the interpretation to be given further on about the 
possible conversion between dark energy and baryonic or dark matter). 
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Regarding the entropy of the universe (equation 10), there are two possible 
interpretations: 

1. Either the entropy of the universe depends only on baryonic and dark 
matter, or the specific entropy of dark energy is zero. In both cases equation 
10 is obtained, therefore, according to the second law of thermodynamics 
S > 0 which implies that (Eq.34) 

TV > 0 (34) 

2. The entropy of the universe is a function of all its components, therefore, 
equation 10 becomes (Eq.35) 

dS = adN + GAdNA (35) 

Where a a is the specific entropy of dark energy and Na is the number of 
dark energy particles. 

Erom the second law of thermodynamics (taking into account that the 
condition Nda + NAdaA = 0 is verified, as deduced from S = Na + Nag a 
and that the universe is an adiabatic system), it follows that (Eq.36) 

S = gN + gaNa > 0 (36) 

Erom equations 32 and 33 different submodels can be deduced, which can 
be classified in two main groups 

5.2.1 T-1 Model 

Particles are created (TV > 0) and dark energy density decreases with time. 
This decrease of dark energy can be interpreted as a consequence of the ex¬ 
pansion of the universe, or as the sum of two processes: the decrease caused 
by the expansion of the universe plus a destruction of dark energy that could 
trigger the creation of baryonic and/or dark matter. 

The entropy of the universe would increase over time if one assumes the 
case of equation 34, or if one assumes the case given by equation 36 as long as 
the condition established by this equation is met. 

By operating in the same way with the Prigogine model, the acceleration 
would be caused according to equation 32 by a negative pressure due to the 
creation of particles. 

A possible interpretation of the evolution of the universe - a qualitative 
one - is shown in figure 1, where the Y-axis represents both terms of equation 
33, and the horizontal axis represents time. In this interpretation there is 
some dark energy in the universe prior to the inflationary and nucleosynthesis 
period. Later, the universe would continue expanding in a decelerated way (all 
densities would decrease - at different speeds - due to the expansion of the 
universe even though there could be particle creation due to the mechanism 
mentioned above). Erom the intersection point of both curves , the universe 
would enter again in an acceleration phase, despite the fact that dark energy 
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Fig. 1 Qualitative representation of both terms of equation 33 (separately, on the Y-axis) 
vs. time (X-axis), arbitrary units in both axes. The almost linear curve represents the dark 
energy component. The intersection point correspond to change in the sign of the accelera¬ 
tion of the universe. 


density would continue to decrease. The subsequent evolution of the universe 
would lead to the elimination of dark energy (in the case where its particles 
are being destroyed), ending the production of baryonic and/or dark matter 
particles. 

5.2.2 T-II Model 

In this model, particles are destroyed {N < 0) and the dark energy density 
increases with time. In this case two situations may arise: 

1. If either the entropy of the universe is only a function of baryonic and dark 
matter, or the specific entropy or dark energy is zero, then this model is 
not possible if the second law of thermodynamics must be fulfilled. The 
other possibility is that this model would indicate that the second law of 
thermodynamics is not true and entropy may decrease with time. 

2. If equation 36 holds, this model could be compatible with the second law 
of thermodynamics. The decrease of baryonic and/or dark matter entropy 
would be compensated by an increase of dark energy entropy. 

In this model, the origin of the acceleration of the universe could be inter¬ 
preted as an anti-gravitational effect or a weakening of the gravitational field 
cause by dark energy (this interpretation would relate to the ideas proposed 
by other authors [27-29]). 

A possible interpretation of the evolution of the universe - a qualitative 
one - is shown in figure 2, where the Y-axis represents both terms of equation 
33, and the horizontal axis represents time. In this interpretation there is also 
some dark energy in the universe prior to the inflationary and nucleosynthesis 
period. Later, the universe would continue to expand in a decelerated way. 
Baryonic energy and dark matter density would decrease (due to the expansion 
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Fig. 2 Qualitative representation of both terms of equation 33 (separately, on the Y-axis) 
vs. time (X-axis), arbitrary units in both axes. The almost linear curve represents the dark 
energy component. The intersection point correspond to change in the sign of the accelera¬ 
tion of the universe. 


and to the elimination of particles) while dark energy density would increase 
(despite the expansion of the universe, this density could increase considering 
that dark energy is the result of the sum of a constant vacuum energy and 
the energy caused by the destruction of particles). From the intersection point 
of both curves, the universe would enter again in an acceleration phase. The 
subsequent evolution of the universe would lead to the complete destruction 
of baryonic and/or dark matter. 


6 Conclusions 

In this research we have obtained a general equation (Eq.l4) that relates the 
variation of the total number of baryonic or dark particles (in a flat, open or 
closed, adiabatic FLRW universe) to the components of the universe (either 
baryonic or dark matter or dark energy). This equation has been obtained 
by balancing equation 8 - equation obtained from the thermodynamic prin¬ 
ciples of adiabatic open systems (and applied to baryonic and dark matter 
and assuming that dark matter will comply with the same laws) - with the 
corresponding part of the equation 13 obtained from the expanded Friedmann 
equations. 

The general condition (Eq.21) that must be met for the universe to show 
stages of accelerated expansion has been deduced, considering the possibility 
of a variation in the total number of particles. 

From the second law of thermodynamics follows that if the entropy of the 
universe is a function only of the baryonic or dark matter, there must be 
particle creation. However, if dark energy may vary over time, then in order 
to fulfill the second law of thermodynamics the condition given by equation 





Title Suppressed Due to Excessive Length 


11 


36 must be satisfied. Also, there is a possibility of particle exchange between 
dark energy and baryonic and/or dark matter. 

From the general equations 14 and 21, two particular cases have been 
studied. The first one is a model originally studied by Prigogine [24] and later 
by other authors [21, 23]. The second case, the Transmutation Model has been 
studied for the first time in this research (as far as the author knows). 

From the Transmutation model follow two possible evolutions of the uni¬ 
verse (in both cases fulfilling the second law of thermodynamics). According 
to the first one (the T-I Model) there would be a creation of baryonic and/or 
dark matter particles from dark energy. The cause of the accelerated process 
would be explained - in a similar way as the Prigogine model - by a negative 
pressure. 

According to the second possibility (T-II model) there would be a destruc¬ 
tion of baryonic and/or dark matter particles resulting in a generation of dark 
energy particles. In this case, the reason for the accelerated expansion of the 
universe should be sought in a weakening of the gravitational force caused by 
the presence of dark energy, or by the possibility that dark energy shows the 
phenomenon of anti-gravitation. 

Taking into account the real lack of knowledge about what dark energy is 
(in case it exists), it is always interesting to deduce new models that predict 
any of its properties. The Transmutation model predicts a particle exchange 
between baryonic and/or dark matter and dark energy, as well as the possi¬ 
ble existence of some anti-gravitational or gravity weakening effect. Particle 
physicists and astrophysicists may have some say in this matter. 
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